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We sliow that tlie key transport states, insulating and conducting, of large regular networks of 
scatterers can be described generically by negative and zero Lyapunov exponents, respectively, of 
Mobius maps that relate the scattering matrix of systems with successive sizes. The conductive phase 
is represented by weakly chaotic attractors that have been linked with anomalous transport and 
ergodicity breaking. Our conclusions, verified for serial as well as parallel stub and ring structures, 
reveal that mesoscopic behavior results from a drastic reduction of degrees of freedom. 
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We describe a remarkable relationship between elec- 
tronic transport in regular arrays of scatterers and a 
special class of low-dimensional nonlinear dynamical sys- 
tems characterized by weak chaos [H Q • This link be- 
tween two disciplines throws light into the nature of 
the insulator-conductor transition in condensed matter 
physics, while the currently studied field of weak chaos 
that exhibits anomalous diffusion and ergodicity break- 
ing is provided with a physical application 21 ■ Quan- 
tum transport properties have become of interest due 
to their fundamental importance in the development 
of nanotechnology; for example, the stability of spin- 
tronic devices based on quantum networks has been in- 
vestigated [3], quantum interference phenomena such as 
Aharonov-Bohm oscillations in conductance, band for- 
mation, and metal-insulator transition in disorderless 
networks, have been studied both experimentally and 
theoretically Quantum networks are also impor- 

tant as theoretical models of molecular devices |8] and 
mesoscopic systems • The latter are experimentally 

available due to the advancement in microfabrication as 
well as to availability of auxiliary tools in the microwave, 
acoustic, elastic, and optical domains Ill4l4|. 



In spite of the importance of understanding the propa- 
gation of electron waves through large networks of quan- 
tum wires with regular or disordered structures, the 
study of the interplay between the individual scatterers 
and their geometric arrangement, that as a whole results 
in complex electronic transport behavior [l5l |. has been 
the object of less exploration. It is, therefore, pertinent 
to develop a procedure to explicitly appraise the inter- 
relations of these two nontrivial facets of the scattering 
processes. A first step in this direction is the determi- 
nation of the generic features of transport due solely to 
the network structural design. Recent advances 15- 17j| 
on the study of coupled limit-cycle and chaotic oscilla- 
tors may serve, unsuspectingly, as a mirror of similar 
simplifying features found here for electronic transport 
properties on large disorderless networks. In the study 
of diffusively-coupled nonlinear oscillators , geometric 



network structures are set up by coupling matrices, and 
the properties of these matrices determine the dynami- 
cal behavior of these systems, independently of the de- 
tails of each individual oscillator. As for globally-coupled 
limit-cycle oscillators, the occurrence of low-dimensional 
nonlinear dynamics in large phase-oscillator systems, ob- 
served sometime ago, has been explained only recently by 
the role that Mobius maps play in controlling the dynam- 
ics of these systems 17]. Since fractional linear transfor- 
mations are invertible, the possibility of chaotic behavior 
would be evidently attributed to the time evolution of 
matrix parameters. 

In this communication we report a general phase tran- 
sition scenario in electronic transport in networks of scat- 
terers connected either in series or in parallel. Specifi- 
cally, we show, by using the scattering matrix approach, 
that the conducting and insulating phases of the meso- 
scopic systems under consideration can be predicted by 
the behavior of the finite-time Lyapunov exponent of a 
nonlinear map in the complex plane, of the Mobius group 
type, that represents the recursive relation between the 
scattering matrices of successive size generations of a net- 
work structure. If we regard the system generation index 
n as the number of map iterations, then the dynamical 
behavior can be used to describe in a quantitative man- 
ner the electronic transport properties of the quantum 
mesoscopic system. 

We consider symmetrical networks constructed by 
putting together a collection of individual scatterers as 
building blocks, shown schematically in Fig. [TJ The as- 
sembled networks have only two end-points and the single 
scatterers when connected in series have two terminals, 
but when connected in parallel have multiple terminals. 
A simple model for three-terminal junctions (connectiv- 
ity = 2) is described by the scattering matrix [l^ 
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where e, a, and /3 are real parameters, related to the 
transmission and reflection amplitudes of the node: < 



2 





FIG. 2. Three different junctions described by an scattering 
matrix St'- (a) a geometric connection, (b) a stub, and (c) an 
Aharonov-Bohm ring threaded by a magnetic field. 



FIG. 1. Construction of regular mesoscopic networks of scat- 
terers in parallel (a) or in series (b) by means of three-terminal 
(a) or two-terminal (b) junctions, where a is the lattice con- 
stant. 



e < 1/2, a = -(l-VT~~2i)/2, and /? = (l + yr~~2i)/2. 
When these nodes are repeatedly connected in paral- 
lel, keeping one initial terminal free, a Cayley tree [2lj 
is formed, and the symmetrical network with two end- 
points consists of a double Cayley tree joined by two- 
terminal individual scatterers described by the matrix 
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which is 2x2, where r?, and 1^ are reflection and transmis- 
sion amplitudes, respectively. Then the scattering ma- 
trix Sn+i of a network of 2(n A- 1) scatterers obtained by 
doubling the size of a previous generation network of 2n 
scatterers with scattering matrix Sn is given by 0] 

-1 



5,, 



n+l 
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where 1 is the 2x2 unit matrix and a is the lattice 
constant. 

Clearly, the network with serial connections is a chain. 
The scattering matrix Sn+i of a chain network of 2(71 -I- 1) 
scatterers obtained by connecting end-to-end two identi- 
cal scatterers Sb with scattering matrix Sn of 2n scatter- 
ers is given by 

^"+1 = . ha - i^l - t^t'b) Sn] ■ (4) 
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For elastic scattering the matrix Sn must be unitary 
as a result of flux conservation, and has the general form 
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which can be diagonalized through the similarity trans- 
formation S'n = USnU\ where U is the unitary matrix 
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and 
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(6) 



(7) 



where On and O'n are the eigenphases, which are related 
to the reflection and transmission amplitudes through 
gje„ = rn + \/tnt'n, e'^" = r„ — yjtnt'n, and consequently, 
the dimensionless conductance (electronic conductance G 
in units of 2e^//i) can be written as g„ — = — 

We rewrite the recursive relations ([3]) and (|H) as one- 
dimensional maps for the eigenphases: dn+i ~ f{Qn) (the 
map for 9'^ is identical) . These maps are fractional linear 
transformations of the form 



Zn+l = F(Zn) = 



AZn + B 



(8) 



e " and A, B, C and D are complex num- 

2ika 



where 2, 

bers: A = 6^^'=°, B = -^1 - 2e, C = -^J\- 2ee2*'=°, 
and D — \ for the double Cayley tree, while for the chain 
network A = -(rg - t'^tb)e^^^°' , B ^ n, C = -ne^^^" , 
and D = \. The transformation F{z) and its inverse are 
analytic on the unit circle in the complex plane, and via 
functional composition define a subgroup (called Mobius 
group) that maps one-to-one the unit circle onto itself. 
The finite n Lyapunov exponent A„ associated to this 
map is defined by 



A„ = - log 

n 
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where 6q is an initial condition. 

The general form ([8]) is due to the symmetry of the 
networks together with the uniform distribution of scat- 
terers, and for this reason, in principle, we can generalize 
the model networks to a class of networks whose growth 
in size is characterized by successive scattering matrices 
generated by Mobius actions. In what follows we report 
numerical results that verify the relationship between the 
conductance of a network and the finite n Lyapunov ex- 
ponent of its associated map f{0n), for which the number 
of iteration time steps is the generation index n that mea- 
sures the size of the network. To demonstrate the gener- 
ality of such a relationship we consider different types of 
junctions. 

(i) A double Cayley tree with connectivity K. The elec- 
tronic transport properties of this structure has been re- 
ported in Ref. for if = 2. Also, for this case, we con- 
sider the following three different junctions (see Fig. [5]) : 
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FIG. 3. The dimensionless conductance g„ (solid lines, n — 
20) and the finite n Lyapunov exponent A„ (dashed lines, n = 
100) are plotted as a function of ka in a network formed by 
a double Cayley tree, when the connecting scattering centers 
are given by (a) a perfect conducting line, (b) a stub, and (c) 
a ring threaded by a magnetic field {L = a, L\ = L2 = L/2, 
$ = 0.1/ic/e). (d) Convergence of the finite n Lyapunov 
exponent A„ to its limiting value A. The sign alternations for 
ka = 1.2 indicate the intermittent quality of the conductive 
phase. 



(a) A geometric connection with the central scatterer de- 
fined by Sb = cr^e^^"', with a Pauli matrix; (b) a stub 
or quantum gate of length L defined by 20] 
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and (c) an Aharonov-Bohm ring threaded by a magnetic 
field whose scattering matrix elements are given by (20l | 
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where Li and L2 are the lengths of the upper and the 
lower arm of the ring, and AL = L2 ^ Li. Here, the 
wave vectors are given by ki = k + 2'iT^e/hcL and 
/c2 = k—2Tr^e/hcL, where L = L1+L2, $ is the magnetic 
flux through the ring, and Afc = k2 — ki = —Ane^/hcL. 

In Fig. [3] we show the relationship between the conduc- 
tance and the finite n Lyapunov exponent for very large 
n. We see that for each case the Lyapunov exponent is 
negative for small values of ka, approaches zero as ka 
increases, and the network is in an insulating phase with 
(7„ = 0. The exponent A„ reaches zero at some critical 
value of ka, and at the same point the mesoscopic system 
undergoes a transition from an insulating to a conduct- 
ing state. As ka is further increased gn > while A„ = 
[\d9n+i/d9o\/n oscillates but becomes zero for n 00, 
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FIG. 4. The dimensionless conductance g„ (solid lines, n = 
20) and the finite n Lyapunov exponent An (dashed lines, 
n = 100) are plotted as a function of ka in a network formed 
by a chain of serially-connected scatterers of (a) stubs, and 
(b) rings threaded by a magnetic field. 



as shown in Fig. ISj^d)], until at a second critical value 
the system returns to the insulating state and A„ < 0. 
There is a remarkable, perfect, correspondence between 
the transport and the dynamical behavior of the nonlin- 
ear maps that represent the three different junctions in 
Fig.Hlvl. 

An extrapolation of the expression for the recursive re- 
lation of scattering matrix, or the eigenphase map f{9n), 
for the general case of arbitrary K, indicates that the 
metal-insulating transitions for K > 2 can be deduced 
from the K = 2 case simply by replacing e by 2e/K. 

(ii) A chain of serially-connected scatterers. The chain 
network when the scatterers are stubs, as in Fig.[2][b), has 
been studied in Ref. 12 by using the transmission matrix 
method. Here we analyze its phase transition proper- 
ties via the scattering matrix approach, and in this way 
reveal some common transport properties of mesoscopic 
systems with differing types of scatterers. Fig. HJa) ex- 
hibits the same transition scenarios as observed in the 
double Cayley tree. A linear network may appear to be 
a more realistic structure than the double Cayley tree, 
but the characteristic behavior the Lyapunov exponent 
appears again to be aprecise indicator of the metal- 
insulating transitions [7'|. For an array of mesoscopic 
rings threaded by magnetic flux as in Fig. [2jc) the el- 
ements of the building block scattering matrix are given 
by Eqs. pT|) . The numerical results shown in Fig. Iljb) 
reveal once more the transition scenarios described pre- 
viously. The positions where the Lyapunov exponent be- 
comes zero indicate quantitatively the locations of the 
metal-insulating transtions that occur in the system of 
serially connected rings. 

Clearly, our results provide convincing evidence that 
low-dimensional dynamical transitions, from a stable 
flxed-point to weakly chaotic attractors, with negative 
and vanishing Lyapunov exponents, respectively, are in- 
dicators of insulating-conducting transitions occurring in 
solid-state model systems. Interestingly, the vanishing of 
the ordinary exponent A takes place via intermittency in 
the sign of the finite n exponent A„ , a signature of weak 
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FIG. 5. The dimensionless conductance Qn (solid) and the 
finite n Lyapunov exponent An (dashed, n = 100) are plotted 
as a function of ka for a chain of delta wells of imaginary 
intensity —iv where v = fi^ /2m. 



apply; the corresponding results are shown in Fig. [5] As 
for asymmetrical arrangements of assemblies of scatterers 
it is expected that chaotic or other complicated scenarios 
may be observed, as the recursive relation might not be 
invertible. 

Furthermore, a significant generic finding, as it is inde- 
pendent of the nature of the network scatterers, is that 
the conducting phase is characterized by the features of 
weak chaos, that is, ergocidity breaking, infinite invari- 
ant measures, and anomalous transport. Experimental 
realizations are feasible nowdays via the use of optical 
lattices as models of condensed matter systems 22| . The 
effect of disorder is left for future studies. 



chaos in maps with a tangency feature like ours 0, 
This constitutes a physical property of the conducting 
networks, i.e. the conductance oscillates when finite size 
is changed. (A unit increment in n represents at least a 
doubling of size in the linear chain). Our development 
also points out a rare simplification in which there is a 
large reduction of degrees of freedom: a system composed 
by many scatterers is described by a low-dimensional 
map. This property is reminiscent of the drastic reduc- 
tion in state variables displayed by large arrays of coupled 
limit-cycle oscillators, for which their macroscopic time 
evolution has been shown 17| to be governed by under- 
lying low-dimensional nonlinear maps of the form of Eq. 
([5]). The parallelism between the transport properties 
of the networks studied here and the dynamical prop- 
erties of arrays of oscillators can be made more specific 
by noticing that in both problems the variables of inter- 
est are comparable, phase shift of the scattered states 
and phase time change of coupled oscillators, both deter- 
mined by the action of the Mobius group. In the scatter- 
ing problem we arrive at the basic nonlinear map by first 
constructing a family of self-similar networks arranged 
by size, or generation n, and then relating the scatter- 
ing matrices of two consecutive generations. In the case 
of coupled oscillators the puzzle of the drastic reduction 
of variables finds a rationale in the identification of the 
role of the Mobius group in the temporal evolution of the 
system. Notice that the particulars of the scattering po- 
tentials appear only through the coefficients in Eq. (jS]) 
while the structure of the network gives the transforma- 
tion its general form. 

Because the basic recursive relation is a consequence 
of special symmetries and an underlying group-theoretic 
structure, it is interesting to explore the relevance of 
our results in other symmetric networks, including two- 
dimensional arrays of interconnected scatterers. We no- 
tice that the correspondence between zero Lyapunov ex- 
ponent and conducting phase does not hold for scattering 
with absorption. In this case the scattering matrix is not 
unitary and, therefore, the Mobius group theory does not 
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